On the puzzle of space weathering alteration of basaltic asteroids 
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ABSTRACT 

The majority of basaltic asteroids are found in the inner main belt, although a few have 
also been observed in the outer main belt and near-Earth space. These asteroids -referred to 
as V-typcs- have surface compositions that resemble that of the 530 km sized asteroid Vesta. 
Besides the compositional similarity, dynamical evidence also links many V-type asteroids to 
Vesta. Moreover, Vesta is one of the few asteroids to have been identified as source of specific 
classes of meteorites, the howardite, eucrite, diogenite achondrites (HEDs). 

Despite the general consensus on the outlined scenario, several questions remain unresolved. In 
particular, it is not clear if the observed spectral diversity among Vesta, V-types and HEDs is 
due to space weathering, as is thought to be the case for S-type asteroids. 

In this paper, SDSS photometry is used to address the question of whether the spectral diversity 
among candidate V-types and HEDs can be explained by space weathering. We show that visible 
spectral slopes of V-types are systematically redder with respect to HEDs, in a similar way to 
what is found for ordinary chondrite meteorites and S-types. On the assumption that space 
weathering is responsible for the slope mismatch, we estimated an upper limit for the reddening 
timescale of about 0.5 Ga. Nevertheless, the observed slope mismatch between HEDs and V- 
types poses several puzzles to understanding its origin. The implication of our findings is also 
discussed in the light of Dawn mission to Vesta. 

Subject headings: Minor planets, asteroids: general - minor planets, asteroids: individual (Vesta) - 
meteorites, meteors, meteoroids 
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1. Introduction 

Vesta is the largest basaltic fully differentiated 
asteroid. Remote visible and near-infrared spec- 
troscopy indicates the presence of basaltic min- 
eralogy on Vesta's surface a nd the possible pres- 



ence of other components jMcCord et al. 1970; 



Larson k Finklll975l : lGaffevl[l997[ ). A large num- 
ber of small asteroids that show a similar surface 
composition, the so-called V- type asteroids , pop- 
ulates the inner m ain belt (jBinzel k Xul 11993 ; 
Duffard et al.l 120041 ) . V-type asteroids have di- 



mensions less than a few tens of km that cannot 
sustain differentiation processes, therefore they 
must have been originated from one -or more- 
much larger, differentiated parent body. Thus, it 
has been su ggested that Vesta is the parent body 
for V-tvoes (jMcCord et allll970h . 
To support this scenario, apart from the aforemen- 
tioned compositional similarity, it has been shown 
that many V- types and Vesta belong to a large as- 



teroid family ( Marzari et al.lll996t INesvornv et al 



2008). To reinforce this link, a 460 km wide crate r 
has been detected on Vesta (jThomas et al.lll997j ). 
Another interesting aspect concerning Vesta and 
V-types is that their visible and near-infrared 
spectra have been found to be similar to those 
of a particular suite of meteorites, namely the 
howard ite, eucrite and diogenite achondrites 
(HEDs:lMcCord et al.lll970l:IConsolmagno k Drake 



19771: lTakedalll997t lDrakell200ll) 

On the other hand, it is also clear that nei- 
ther all V-types nor HEDs are compatible with 
being derived from Vesta; therefore, it is sus- 
pected that they origina te from multiple differ 



entiated parent bod ies (jYamaguchi et al.l 12002 ; 
Wiechert et al.l [2004 ) . In this respect, it is sig- 



nificant that a few V-types have been discov- 
ered beyond the 3:1 mean motion res onance 
(jLazzaro et al.ll2000t iDuffard k Roigll2009h which 
cannot be related to Vesta on dynam ical basis 
(|Roig et all 120081: INesvornv etaH I2OO8I ) and -at 
least i n the case of Magnya- also on compositional 
basis (jHardersen et al.l 120041 ) . 
Despite the general scenario described above is 
largely accepted, there are some observations 
that remain unexplained. First of all, V-types 
generally exhibit stronger absorption bands than 
Vesta and steeper con tinuum both in the visible 
(Ifflroi k Pieterslll998l) an d in the 0.4-1.6 fim re- 
gion ([Burbine et al.ll200lh . 



The reason for these discrepancies is not yet 
clear. Possible causes are indicated in compo- 
sition/texture differences or due to the effects o f 
some alteration process (e.g. Burbine et al.ll2001l) , 
like for instance those that have been demon- 
strated to cause optical alteration on silicate-rich 
S-type asteroids (|Jedicke et al . 2004; iMarchi et al 
2006a|). However, on S- types the 1 and 2 ^m ab- 



sorption bands tend to become shallower by in- 
creasing space weathering (i.e. increasing spectral 
slope), while this effect is not observed on V- 
types. V-type asteroids that have reddest slopes 
also show strong 1 an d 2 /im absorption bands (e.g. 



Burbine et al.l 120011 ). Vesta is found to be most 



comparable to h owardites with a fin e-grained size 
distribution (e.g. Pieters et al.ll2006h ■ 
On the other hand, the laboratory experiments 
conducted on HEDs and pyroxenes confirm that 
V-type materials should alter under the effects 
of ion bombardment and micrometeorite im- 
pacts, in a way similar to S-types, although 
the al teration timescale is likely somehow re- 
duced dHiroi fc Pietersl[l998l: IMarchi et al.ll2005al 



Vernazza et al. 



20061 ). 



In this paper we present a comprehensive model 
of space weathering of candidate V-types detecte d 
via SDSS photometry (|Roig k Gil-Huttonl [2006^ . 
The present analysis aims to investigate whether 
or not V-type surfaces are affected by space weath- 
ering processes. We identify several interesting 
facts that support the presence of space weath- 
ering alteration, although its behavior appears 
distinct to what found on S-types. 
In doing so, our findings disclose interesting as- 
pects of surface properties of Vesta and V-types 
that are useful also in the context of the NASA 



Dawn mission (jRussell et al.l I2007T ) for what con- 
cerns mission planning and Vesta's data interpre- 
tation. 

2. Modeling the space weathering of V- 
type asteroids 

Spectrophotometric catalog of V-types 

It has been demonstrated for S-type asteroids that 
visible spectral slopes can be successfully used to 
model the optical alterations of asteroidal surfaces 
due to space weathering. Spectral slopes have 
been evaluated using di fferent methods, m ainly 
photometric data (e.g. Jedicke et al. 20041 ) and 
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spect roscopy (jMarchi et al.ll2006ai iPaolicchi et al 
20071 ). leading basically to the same conclusions. 
Concerning V- types, spectroscopic data are avail- 
able only for a restricted number (~70) of ob- 
jects, preventing any significant characteriza- 
tion of the whole population. For this rea- 
son, in the present work, we will make use of 
SDSS photometry reported i n the Moving Ob- 
ject Catalog V3 (MOCV3; llvezic et all 12001 



Juric et all l2002h . which allows the analysis to 
be extended to include a large number of ob- 
jects (388 and 466 family and non-family V- 
types, respec tively). V- types ha v e bee n selected 
according to Roig fe Gil-Hutton (2006). In ad- 
dition, we also include 52 V-types identified by 
the spec troscopic surveys SMASS, S3QS2 and 
SINEO dLazzaro et alj|2004 bus fe Binzell 12001 
Binzel et alJl2004t iMarchi et al.ll2005bl ) which are 
not included in the SDSS MOCV3. It is im- 
portant to underlin e that the list established by 
Roig fe Gil-Huttonl ()2006h consists of candidate 
V-types and that spectroscopic data are needed 
for a robust identification. Nevertheless, a com- 
parison between candidate V-t ypes and confirmed 
ones displays a good match ( Roig fe Gil-Huttonl 
20061) . 

Spectral slopes have been evaluated by linear best 
fit of the albedo using the SDSS bands g', r',i', 
z'. One sigma errors of SDSS albedoes have been 
used as weighting factors for the fitting procedure. 
As for spectroscopic data, we first evaluated the 
albedo in g', r', i', z' bands through cubic spline 
of the spectra and then proceeded with linear best 
fit. Albedoes have been n ormalized to the r' band 
(iRoig fc Gil-Huttonl 120061 ). 



Parameterizing the space weathering 

Previous analysis of S-types and ordinary chron- 
drites (OCs) showed that the space weathering 
on asteroids is a process that, if no other alter- 
ation is present, is expected to progressively evolve 
over time until eventually a sort of saturation is 
reached. Moreover, its efficiency depends on the 
location and past orbital evolution of the aster- 
oids. 

In this work, the study of space weathering is 
done according to the model dev e loped in a se- 
ries of papers by Marchi et al. ( 2006al fbT) and 



that the aging, location, and past evolution can be 
condensed into a single parameter, named expo- 
sure E. The exposure is proportional to the dose 
of radiation/ion received from the Sun and scales 
as 



E = 



T 



2 VT 



(1) 



where a, e are the Keplerian eleme nts of the 
body, and T is its average collisional age (jMarchi et al 
2006a)). For MB As we used proper elements given 
their stability on a longer timescale. Ages are de- 
rived according to the collisional evo l ution of main 



belt asteroids (Bottke et al. 2005; Marchi et al 



2006aD . Collisional lifetimes depend upon the di- 
ameters, which have been estimated using the 
measured avera ge geometric albedo ( 0.3) of a sam- 
ple of V-types ( Tedesco et al. 20021 : iBenner et al 



20021: iDelbo et al.ll2003ll2006l) . The collisional ag. 



of V-types spans from about 0.5 to 2.7 Ga. These 
extremes correspond to the smallest and largest 
bodies of the sample, which have estimated diam- 
eters of about 1 km and 10 km, respectively. In 
what follows, it is assumed that the formation of 
Vesta's family occurred sometime before the old- 
est age in the sample, namely ~ 2.7 Ga, so that 
asteroid's age are set by subsequent collisional 
evolution. Although Vesta's fa mily age is not 



well c onstrained, a recent work of iNesvornv et al 



(2008) bund evidence for a much earlier origin, 
possibly dating 3.5-3.8 Ga a go, than previously 
assumed ( Marzari et al. 19961 ). 
For the NEOs, given their short evolution timescale 
in near-Earth space (~ 5 Myr), the major contri- 
bution to their space weathering alteration comes 
from their past evolution into the main belt. 
Therefore, for each NEO we set a equal to its 
compu ted average source distance ()Marchi et al 
2006al ) and e = 0.17 (the average of MBAs). 



Paolicchi et al.l ([20071 ). This model demonstrated 



3. Space weathering on V-types? 

The main results of our model can be inves- 
tigated with the help of the slope-exposure scat- 
ter plot. For this purpose, we separated V-types 
belonging to the Ves ta dynamical family from 
non- family members ( Roig fc Gil-Hutton 2006). 
Results are shown in Figure [TJ For comparison, 
the trend derived for S-types is also shown. Note 
that the latter is computed by the least-square 



3 



fitting of S-type visible spectra of MBAs detected 
by SMASS and S30S2 surveys (a total of 750 ob- 
jects), using the SDSS bands fitting procedure, 
as described above. This trend is very close to 
th e one derived by mea ns of spectroscopic data 
in IPaolicchi et all (|2007l) (valid for E > 10), en- 



suring that the procedure adopted here is robust. 
Figure Q] also reports V-type NEOs along with av- 
erage slopes and range of slopes for a sample of 50 
HEDs. Note that these HED samples have similar 
grain size (< 25 /im); thus, the observed spread 
in slope is likely due to differences in composition, 
rather than due to grain size effects. 
A number of interesting results can be derived 
from Figure [Tj First of all, V-types are, on aver- 
age, redder than HEDs. This difference is mini- 
mized for eucrite meteorites. To better study this 
difference, we report the slope histograms in Fig- 
ure [21 The average slope gap between eucrites and 
MBA V-types is abo ut 0.6 Mm" 1 , similar t o that 
of S-types and OCs (|Paolicchi et al.ll2007l) . The 



largest difference, of about 1.0 /im -1 , is found for 
diogenites and MBA V-types. 
It is not easy to understand the origin for this 
marked difference with available data. It could 
be due, for instance, to systematic differences in 
grain size or observational conditions, such as the 
phase angle. However, it is hard to believe that 
these factors can be entirely responsible for the ob- 
served slope variations for the following reasons. 
Concerning the grain size, we computed the effect 
of the particle size on the slope by analyzing two 
eucrites and one diogenite for which spectra with 
different particle sizes (in the range 25-250 /zm) 
have been gathered (RELAB database). In this 
range, the slope increases with de creasing particle 
sizes by about 0.4 Mm -1 ( see also lHiroi fc Pietera 



1998; iBurbine et all 120011 ). Therefore, to obtain 



the observed reddening on asteroids, their parti- 
cle size should be much smaller than 25 /im, a 
conclusion which seems unrealistic given the low- 
gravity environment of V-typetQ. Moreover, we 
recall that the HED slopes used here are virtually 
the reddest, given that the particle size of these 
samples is below 25 /im. Concerning the phase 



1 For comparison, the average regolith particle size on the 
Moon is 70 /im, and only 10-20% of the particles have sizes 
< 20 fim. Moreover, thermal inertia modeling shows that 
asteroid's regolith is m uch coarser than the lunar one (e.g. 

iDelbo' fc Tangall2009l 1. 



angle of observations, it spans from about 2 to 
30 deg, which is not too different from geome- 
try condi tions during spectra acquisition of HED 
samples (jBurbine et al.l l200lh . There also seem 
to be differences between family and non-family 
V-types. The latter shows a larger slope scatter 
and has a considerable number of objects having 
slope < —0.6 /im -1 and > 0.5 /im -1 (see Fig. [2]). 
These two extreme classes of object are separated 
by the bulk of the slope distribution, and might 
represent two distinct populations of objects. The 
"blue" non- family V-types (< —0.6 /im -1 ) are 
compatible with HEDs, and do not show any ef- 
fect of spectral alteration. 

The "red" V-types (> 0.5 /im -1 ), are much redder 
than HEDs. Their exposure is not systematically 
different than other V-types; therefore the ob- 
served red slope is not justifiable in terms of space 
weathering. Possible interpretations are that: ei- 
ther i) they are S-types that are erroneously clas- 
sified as V-types or ii) if basaltic, they have a pe- 
culiar composition which is not sampled by HEDs. 
After a close look at SDSS data for these objects, 
it seems that i) could be ruled out; however, v-nir 
spectroscopy is needed for a secure classification 
(see end of this section for further comments on 
this issue). Moreover, these extreme slope objects 
do not have any significant difference in terms of 
orbital parameters with respect to other V-types. 
Although less evident, blue and red objects are 
also present among family V-types. 
Interestingly, NEO V-types are somehow in be- 
tween HEDs and MBAs. This is, at least for 
some NEOs, in agreement with their computed 
i ntermediate exposure , as also found for S-types 
([Paolicchi et al.l 120071) . Nevertheless, the optical 
properties of NEOs can be a ffected by other pro- 
cesses like close e ncounters ( Marchi et al. 2006b; 
Binzel et alll2010l ). The effects of the latter on V- 



type NEOs will be investigated in a future work. 
Alternatively, several differences with S-types can 
be noted. The most striking is that the slope 
shows a negative trend for increasing exposures. 
In other words, the more exposed objects tend 
to appear bluer, on average, than less exposed 
ones. The trend is less pronounced for non-family 
members, and becomes remarkably noticeable for 
family members. A linear best fit performed 
through minimizing chi-squared error statistics 
gives the following trends (expressed in units of 
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slope[/xm- 1 ]/log(£;[MaAU- 2 ])): -0.71 ±0.02 and 
—0.35 ± 0.02 for family and non-family V-types, 
respectively. Error bars are estimated taking 
into account 1-sigma errors of slope values. Both 
trends are statistically highly significant with a t- 
distribution two-tailed probability of significance 
of 10~ 6 and 10~ 9 , respectiveljU. For comparison, 
the S-types trend is 0.23 ± 0.07 (two-tailed prob- 
ability of 10~ 3 ). 

In light of the above results, the unexplained lack 
of spectral reddening on Vesta can be studied 
from a new perspective. First of all, there seems 
to be a general de-reddening trend with increas- 
ing exposures. These trends are also marginally 
present in terms of diameters, and become pro- 
gressively more pronounced in terms of ages and 
exposures. This result suggests that the trends 
are not merely due to an asteroid size effect. Note 
that Vesta being merely the largest body (i.e. the 
most exposed) fits with this trend. 
The slope scatter of V-types is larger than that 
of HEDs, however the latter may be limited due 
to the relatively low number of samples used in 
this work. Therefore, as for the HEDs, it is pos- 
sible that much of the V-type slope spread may 
be due to compositional/texture differences. It 
could also be (partially) due to the presence of 
S-type -erroneously classified as V-types- interlop- 
ers within the lists used in this work. S-types are 
expected to have, on averag e, higher slopes than 
V-types (|Bus fc Binzelll2002l ). Therefore, the pos- 
sible presence of S-type interlopers may affect the 
observed excess of red objects, which likely has a 
negligible influence on the clear slope gap between 
HED and candidate V-types shown in Fig. [2j 



4. Discussion and Conclusion 

The main result of this paper is that V-type 
slopes are systematically redder than HEDs (see 
Fig. [U Fig. [2]). Thanks to a large sample of can- 
didate V-types, the present work considerably ex- 
tends previous resu lts (e.g. Hiroi fc Pietersl [l998t 
Burbine et al.l 120011 ) and allows quantitative esti- 
mates of the reddening of V-types. The spectral 
mismatch mimics what was found for OCs and 



crepancy could be due to factors such as system- 
atic compositional/texture differences, different 
conditions of observations or space weathering. 
We showed that observational conditions and tex- 
ture are unlikely to be the source of the observed 
discrepancy. Therefore, either the marked slope 
shift between HEDs and V-types is due to sys- 
tematical differences in composition, or it can be 
ascribed to some kind of space weathering alter- 
ation. The previous hypothesis would imply that 
the suite of HED compositions is not representa- 
tive of the whole V-type population. While this 
is certainly an open possibility, it is mandatory 
to obtain v-nir spectra of small V-types to con- 
firm. In the case that the slope mismatch is due to 
space weathering processes, it is possible to con- 
strain the reddening timescale on the basis of the 
collisional ages. Although a precise estimate is not 
possible with present data, an upper limit for the 
reddening timescale is set by the youngest V-type 
MBAs which have ages of ~0.5 Ga. Moreover, 
58% of V-type NEOs have a slope compatible 
with HEDs (within 1-sigma of average values), 
while only about 20% of S-type NE Os are com- 
patible with OCs (|Binzel et al.ll2010h . Therefore, 
if no other effect comes into play, it can be in- 
ferred that V-type's space weathering timescale is 
larger than a factor of ^3 with respect to that of 
S-types. 

The possibility that space weathering is respon- 
sible for the spectral mismatch is also reinforced 
by the few experiments published to date of ion 
bombardm ent of Bereba eucrite and Johnstown 
dioge nite ( Hiroi fc Pietersl 1998 : IVernazza et al 



S-types (jPaolicchi et al.l 120071 ). The observed dis 



2 A two-tailed probability less than 10 2 indicates full sta- 
tistical significance of the correlation. 



20061 ). They both exhibit reddening (and darken- 
ing). Note that these experiments also confirm a 
slower timescale with respect to OCs. 
The true nature of the spectral mismatch still re- 
mains elusive in several aspects; nevertheless some 
noteworthy characteristics have been identified. In 
particular, we explore in detail the possibility that 
the spectral mismatch is due to space weathering 
processes. 

First of all, the spectral slope anti-correlates with 
exposure to solar wind. This is the opposite of 
what was found for S-type s and other spectral 
types ( Lazzarin et al. 2006[) . Therefore, if some 
kind of Sun-related space weathering is operating 
on V-types, either it behaves in a different manner 
with respect to S-types or other factors wipe out 
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the reddening-exposure relation. Concerning the 
latter, a possibility is represented by a composi- 
tion gradient across V-types, from eucrite-like for 
the smallest members (~1 km) to diogenite-like 
for the largest members (~10 km). This puta- 
tive composition gradient would help to explain 
the observed slope-exposure trend for two reasons: 
i) diogenites are bluer, on average, than eucrites; 
and ii) according to the few experiments available, 
the reddening of diogenites seem less pronounced 
than for eucrites, although the space parameters 
of HED compositions and alteration processes are 
far from being exhaustively investigated. 
Note that the above explanation is preferred with 
respect to other possibilities, e.g. grain size vari- 
ation, since the latter would require a grain size 
finer than ~ 25 /im (i.e. higher slope) for smaller 
objects. This is the opposite of what is expected, 
since smaller objects tend to have a larger grain 
size due to the low gravity. Also, the substantial 
presence of S-type interlopers among the red V- 
types would not greatly affect the observed slope- 
exposure trend since they are nearly equally dis- 
tributed in terms of E (see Fig. Q]). 
If confirmed, the above compositional trend 
among family V-types could be the result of cra- 
tering ejecta from various depths in Vesta's crust 
(eucrite is expected in the upper crust, while dio- 
genite is expected in the lower crust); or may be 
due to different episodes of cratcring occurring 
in regions with different compositions. However, 
since large fragments are expected to originate 
in the near-surface spall region, the second hy- 
pothesis seems more probable. This may also be 
confirm ed by the variety 39 Ar- 40 A r shock ages of 
HEDs (|Bogard k Garrison! 120031 ) . Furthermore, 
the spectral data of V-types see m to indicate 



that both mineralogies are present (D uffard et al 



20041 ). but the limited sample prevents a firm con- 
clusion. 

The proposed scenario could be tested as soon as 
more v-nir spectra of small V-types would become 
available. Note that the above discussion holds 
for both family and non-family V-types. For the 
latter, however, the anti-correlation is less pro- 
nounced, and it becomes nearly flat if the outliers 
(red and blue) are excluded. This could be an 
indication that most non-family V-types have a 
different origin with respect to family V-types. 
The NASA Dawn mission will aid in understand- 



ing the weathering process on Vesta, thanks to 
spectral imaging and global mapping. Therefore, 
the possibility to have detailed spectral informa- 
tion of specific areas of different composition and 
texture will shed light on the different effects of 
space weathering on basaltic asteroids. 

We thank the anonymous referee for the helpful 
comments on the manuscript. 
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Fig. 1. — Slope vs exposure scatter plot for MBA V-types. Both family and non- family V-types are reported. 
The position of Vesta is indicated by the plus sign (upper panel) . Horizontal lines indicate the average values 
for a sample of 14 howardites, 26 eucrites, 10 diogenites and 180 OCs. The range of slope values for HEDs 
(computed as 1-sigma of the slopes) is also reported (placed arbitrarily at E = 8,9,10). NEOs slopes are 
indicated by filled green circles. Red dashed line indicates the reddening trend for S-types, while blue lines 
indicate MBA V-types best fit. The average 1-sigma slope error is reported at the upper-right corner of both 
panels. 
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Fig. 2. — Histograms of visible slopes for family and non-family V-types. Average slope values of HEDs are 
also indicated by vertical lines (see also Figure Q}. Note that both distributions show marked non-Gaussian 
shapes, with a clear excess of "blue" and "red" objects, i.e. having slopes < —0.6 /im" 1 and > 0.5 /im _1 , 
respectively (see the text for details). 
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